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We study the viscous damping of r-modes of compact stars and analyze in detail the regions 
where small amplitude modes are unstable to the emission of gravitational radiation. We present 
general expressions for the viscous damping times for arbitrary forms of interacting dense matter 
and derive general semi-analytic results for the boundary of the instability region. These results 
show that many aspects, such as the physically important minima of the instability boundary, are 
surprisingly insensitive to detailed microscopic properties of the considered form of matter. Our 
general expressions are applied to the cases of hadronic stars, strange stars, and hybrid stars, and 
we focus on equations of state that are compatible with the recent measurement of a heavy compact 
star. We find that hybrid stars with a sufficiently small core can "masquerade" as neutron stars and 
feature an instability region that is indistinguishable from that of a neutron star, whereas neutron 
stars with a core density high enough to allow direct Urea reactions feature a notch on the right 
side of the instability region. 



I. INTRODUCTION 

Pulsars rotate with periodicities whose stability ex- 
ceeds that of any terrestrial clock. Pulsar frequencies and 
their time derivatives are thereby by far the most accu- 
rately measured properties of compact stars, whereas all 
other information about them is subject to much greater 
uncertainties. Correspondingly it is extremely tempting 
to exploit this information in order to learn about their 
internal structure and in particular whether they consist 
of novel phases of dense matter that might contain decon- 
fined quarks [T] . This requires the development of unique 
signatures that connect particular microscopic properties 
to the macroscopic data. 

Pulsar frequencies can change over time both by ac- 
cretion of matter that transfers angular momentum from 
a companion star and by the emission of gravitational 
radiation due to oscillations of the star. A particularly 
interesting class of oscillation modes are r-modes [2l-|4] 
which are counter-rotating modes of a rotating star and 
are in the absence of viscous damping unstable at all 
rotation frequencies jS]. This instability transforms ro- 
tational into gravitational wave energy and leads to an 
exponential rise of the r-mode amplitude. When viscous 
damping is taken into account the star is stable at low fre- 
quencies but there remains an instability region at high 
frequencies jFl H] • If this instability is stopped at a large 
amplitude, r-modes are a strong and continuous source 
of gravitational waves and could provide an extremely ef- 
ficient mechanism for the spin-down of a young compact 
star Observational data for spin frequencies of pul- 

sars, that spin down and allow the determination of an 
approximate age associated to their spin-down rate^, is 



The spin-down age is determined by the assumption that mag- 
netic dipole breaking, which features a qualitatively similar be- 
havior as gravitational wave emission, dominates the spin-down. 
It can thereby only give a rough order of magnitude estimate 
for the age. Age estimates are, in particular, not available for 



shown in fig. [T] Whereas observed old pulsars in binary 
systems can spin nearly as fast as the maximum Kepler 
frequency, above which the binding force cannot counter- 
act the centrifugal pseudo-force anymore, and can feature 
rotation periods in the milli-second range, younger stars 
are far below this limit. This is surprising since in their 
creation during a supernova a significant fraction of the 
angular momentum of the initial star should be taken 
over by the much smaller compact core which therefore 
should dramatically spin up. This naive assumption is 
backed up by explicit analyses where millisecond rotation 
frequencies at birth are indeed possible [1^ . If r-modes 
spin down compact stars on time scales shorter than the 
age of the youngest observed pulsar, the lower bound- 
ary of the instability window should give an upper limit 
for the maximum rotation frequency of young compact 
stars^. 

R-modes are also relevant for the case of older stars in 
binaries that are spun up by accretion since they gen- 
erally limit the maximum possible rotation frequency 
of a star to values substantially below the Kepler fre- 
quency. A challenging finding is that, in contrast to 
purely hadronic stars, more exotic possibilities like self- 
bound strange stars [TT], hybrid stars [7] or stars where 
hyperons are present in the core [H] can feature so-called 
"stability windows" where over a range of intermediate 
temperatures the r-mode instability is absent up to rather 
high frequencies. The observation of stars rotating at 
such frequencies could therefore provide evidence for ex- 
otic phases in their interior. In this context the masses 



stars that currently spin up and which are correspondingly not 
included in the plot. However, for the youngest stars there are in- 
dependent age determinations from the observation of the corre- 
sponding supernova remnant that qualitatively agree with these 
estimates. 

^ Note that the final frequency of the spin-down evolution can lie 
below the minimum of the instability boundary since the damp- 
ing of the r-mode can take some time to complete, even after it 
enters the stable regime. 
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and radii of stars provide further important information. 
The recent precise measurement of a heavy compact star 
with M « 2Mq |13I puts constraints on the presence 
of exotic phases since such phases lead to a softening of 
the equation of state which in general leads to a smaller 
maximum mass that is achievable for such an equation of 
state. In combination with pulsar data this should lead 
to more restrictive bounds on the possible presence of 
certain forms of matter in compact stars. 

A major problem for the extraction of information on 
the composition of compact stars from observational data 
is the huge theoretical uncertainty in the equation of state 
of dense matter and its transport properties. This holds 
both for the hadronic side, where nuclear data is only 
available at low densities and large proton fractions, and 
also for hypothetical phases of quark matter, since QCD 
as the fundamental theory of strong interactions cannot 
be solved so far in this non-perturbative regime. An- 
other big uncertainty factor is the crust of the neutron 
star since although its microscopic physics is in general 
more constrained by experimental data than the high 
density phases in the core, its structural complexity lim- 
its so far a complete description. Simplified estimates 
suggest that the crust could have a strong impact on 
the damping of r-mode oscillations via surface rubbing at 
the crust-core interface [15]. However, as in all present 
analyses, r- modes are considered as solutions to the hy- 
drodynamics equations of an ideal fluid. When certain 
regions feature viscosities that are enhanced by orders of 
magnitude, like at the crust-core interface, the r-mode 
profile should strongly change in these regions and the 
reduced amplitude there would result in a considerably 
weaker damping. Therefore, we neglect crust effects in 
our present analysis, but a more complete understand- 
ing of these effects in the future is desirable. Unfortu- 
nately due to all this, even if two phases feature signifi- 
cant qualitative differences these are often overshadowed 
by the huge quantitative uncertainties in the detailed mi- 
croscopic properties of either of them. However, it was 
previously observed that certain features, like the impor- 
tant case of the minimum of the instability region can 
be surprisingly insensitive to quantitative details of the 
considered models [6 . If such statements can be sub- 
stantiated this could allow us to devise robust signatures 
of the qualitative features that can be stringently tested 
with present and forthcoming astrophysical data. 

To this end, we study in this paper the instability re- 
gions of small amplitude r-modes in detail. In contrast 
to previous treatments that studied particular star mod- 
els numerically, we derive general analytic results that 
are valid for stars consisting of various forms of matter 
and compare these to numeric evaluations for realistic 
equations of state. We find that although the form of 
the instability region can be qualitatively distinct for the 
different forms of matter, many aspects of the instabil- 
ity regions are extremely insensitive to the detailed un- 
known microscopic input, like the transport coefficients 
of a given phase. Moreover, we reveal the parametric 
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Figure 1: Rotation frequencies of observed pulsars versus their 
approximate spin-down age from the ATNF pulsar catalogue 



dependence on the underlying microscopic parameters 
where such dependences are significant. We study ex- 
plicitly the cases of neutron stars, strange stars as well 
as hybrid stars and in view of the recent discovery of a 
2 Mq compact star, we generally apply only equations of 
state that can accommodate such a heavy star. We find 
that the instability regions of the hybrid stars are almost 
indistinguishable from those of neutron stars if the size 
of the quark matter core is smaller than roughly half of 
the star's radius. Further, we also study the case of an 
ultra-heavy neutron star where direct Urea reactions are 
allowed in the core and find that it features a notch at 
the right hand side of the instability region. Finally, in 
addition to the dominant m = 2 r-mode we also con- 
sider higher multipoles and note that they could easily 
be excited and become important for the evolution of the 
star. 



In a companion paper |16| we show that due to the 
strong increase of the bulk viscosity with amplitude [17|- 
20 the r-mode instability is only present at sufficiently 
small amplitudes and the exponential r-mode growth is 
eventually saturated at finite amplitudes so that r-modes 
could indeed provide a viable mechanism for both the 
spin-down of young stars and the frequency limit of old, 
accreting stars. 



3 



II. STAR MODELS AND R-MODES 

A. Static star models 

The analysis of compact star oscillations and their 
damping requires as a first step the stable equilibrium 
configuration of the star. In this section we will dis- 
cuss the considered star models that are used later on. 
The equilibrium star configuration is determined by the 
gravitational equations for a fluid sphere. We employ 
here the general relativistic Tolman-Oppenheimer-Volkov 
(TOV) equations |22] . The latter require the equation of 
state of neutral and /3-equilibrated dense matter. The re- 
cent measurement of a neutron star with the large mass 
1.97 ± 0.04 Mq [HI E] puts bounds on the equation of 
state of dense matter. We generally study equations of 
state that can accommodate such a heavy star. We con- 
sider three qualitatively different classes and study in 
each case both a star model with a standard value of 
1.4 Mq and one with a large mass 2 Mq. 

1. Neutron stars (NS) are obtained for an equation 
of state that is hadronic at all relevant densities. 
Whereas the above maximum mass does not pose 
problems for neutron stars obtained from most 
hadronic equations of state, it seem nearly impos- 
sible to obtain such heavy stars which contain a 
significant amount of hyperonic matter [T3] . There- 
fore we do not study this possibility here and con- 
sider only stars consisting of neutrons, protons 
and electrons as well as muons at sufficiently high 
density. We also neglect the possible presence of 
hadronic pairing which in general significantly re- 
duces the bulk viscosity [23', 24J. Such pairing is 
only realized in certain shells within the star and 
as long as they are not large this might not qual- 
itatively change the results obtained here. How- 
ever, the intricate dynamics of a two-component 
fluid can change this simplistic picture [55] . For our 
numerical results we employ the equilibrated equa- 
tion of state by Akmal, Pandharipande and Raven- 
hall (APR) PS] which relies on a potential model 
that reproduces scattering data in vacuum supple- 
mented by a model for three-body interactions in 
order to reproduce the saturation properties at nu- 
clear densities. As a low density extension of the 
APR data we use P71 128j . Furthermore, we study 
in this class also a neutron star with an ultra-high 
mass 2.21 Af0 close to the mass limit for the APR 
equation of state, since in this case direct Urea in- 
teractions are possible in the interior of the star 
leading to a signiflcantly enhanced weak rate, cf. 

e.g. m- 

2. Hybrid stars (HS) with an outer hadronic part and 
a core of quark matter are obtained from an equa- 
tion of state where at some density the effective 
degrees of freedom change from hadrons to quarks. 
In general the equation of state of interacting quark 



matter is unknown and there are only hints from 
the perturbative regime t29ti31j or model studies. 
We use the simple quartic parameterization for the 
equation of state of ungapped 3-flavor quark matter 
[17| [32j in terms of the individual quark chemical 
potentials /id, /i„ and //^ and the electron chemical 
potential /^e 
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where c, nis and B are effective model parameters 
that incorporate some effects of the strong interac- 
tions between the quarks. From this equation of 
state we flnd the /3-equilibrated and charge neutral 
ground state which depends on a single quark num- 
ber chemical potential fig . We use the general form 
since the computation of transport properties be- 
low requires susceptibilities around the equilibrium 
state. Within the parameterization eq. ([T|) the re- 
cent measurement of a heavy star strongly restricts 
the equation of state so that only equations of state 
that are strongly interacting (c > 0.3) are compati- 
ble [m |32] . These are equations of state where the 
transition to quark matter is at rather low values 
of the baryon number < 1.5 no, where uq is nuclear 
matter density, and we choose here one where it 
occurs at 1.5 uq. Since the APR equation of state 
happens to be very similar to the above form eq. 
dl| so that even multiple transitions are possible 
|32| , we do not expect the transition density to be a 
robust result that is independent of the considered 
equations of state. Therefore we study here in ad- 
dition also two 1.4 Mq hybrid star models obtained 
with quark equations of state that cannot accom- 
modate a heavy star when combined with the APR 
equation of state. Ona has a small quark core, ob- 
tained for a transition density of 3.25 tiq, and the 
other a medium-sized quark core, obtained for a 
transition density of Suq. 

There are many possible phases of quark matter 
that feature various color superconducting pairing 
patterns. Here we do not explicitly study star mod- 
els with pairing, but note that the parameterization 
eq. ([T| can also describe superconducting matter. 
For our hybrid star models we make the assump- 
tion of local charge neutrality, excluding the possi- 
bility of a mixed phase and thereby circumventing 
the description of the wealth of possible geometric 
structures of such a mixture. 

3. Strange stars (SS) that are self-bound could ex- 
ist according to the strange matter hypothesis [33] 
that the true ground state of strongly interacting 
matter is 3-flavor quark matter. For strange stars 
we use the same equation of state eq. ([ij, but 
use parameter sets that realize the strange mat- 
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ter hypothesis. In contrast to the case of hybrid 
stars large mass strange stars are possible even for 
c = and due to our ignorance of the precise 
form of the interacting equation of state we choose 
this value to keep our model as simple as possi- 
ble. For a strange quark mass of iris = 150 MeV 
stable strange stars exist in this case for bag con- 
stants Bum < (158 MeV^)'', whereas heavy strange 
stars exist for this mass value only for bag constants 
B < (UOAIeV)'^. Lower effective quark mass val- 
ues or changes in the quartic term (c ^ 0) relax 
these bounds. 

We note already at this point that although the detailed 
parameters we chose here for our star models are rather 
arbitrary, we will give general analytic expressions below 
that will reveal the dependence of our main results on the 
various model parameters. The characteristic parameters 
of the considered star models are given in table |l] As is 
a well known property, the radii of the different models 
vary only moderately with the mass for masses of 1 to 
2Mq. The density profiles are shown in fig. [2] The core 
densities of these star models range from below 3 no to 
more than 7 no- In contrast to neutron and hybrid stars 
that vary in density by 14 orders of magnitude, strange 
stars feature a roughly constant density profile. 





M[Mq] 


M^ore [Mq] 
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Qk [kHz] 


NS 


1.4 


(1.39) 


11.5 


3.43 


1.58 


6.02 
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(1.99) 


11.0 


4.91 


2.46 


7.68 




2.21 


0.85 


10.0 


7.17 


3.37 


9.31 


SS 


1.4 




11.3 


2.62 


1.91 


6.17 




2.0 




11.6 


4.95 


2.43 


7.09 


HS 


1.4 (S) 


0.38* 


10.8 


5.89 


1.85 


6.61 




1A{M) 


0.66* 


10.3 


6.66 


2.09 


7.06 




1.4(L) 


1.06 


12.7 


2.32 


1.17 


5.16 




2.0 


1.81 


12.2 


4.89 


1.84 


6.62 



Table I: Results of the considered models of neutron stars 
(NS), strange stars (SS) and hybrid stars (HS). Shown are 
the mass of the star M, the mass of the core A'Icore, the ra- 
dius R, the baryon density at the center of the star ric given 
in units of nuclear saturation density no, the average density 
(n) and the Kepler frequency Qk- The neutron stars were ob- 
tained by solving the relativistic TOV equations for catalyzed 
neutron matter using the APR equation of state |26| with low 
density extension [271 128) and the strange stars with a quark 
gas bag model with c = 0, nia — 150 AdeV and a bag pa- 
rameter B = (138 MeV)*. Large mass hybrid stars are only 
found when strong interaction corrections are considered, cf. 
[5^ . and we find a 2M0 star for c = 0.4, rUs = UOAleV, 
B = (137 MeV)". The additional two 1.4 Mq hybrid mod- 
els with smaller cores, marked with an asterisk, result from 
equations of state that do not allow large mass models. They 
correspond to c = 0, = 150 A/eV, B = (164.5 MeV)" and 
(171.5 MeV)* which are chosen to obtain transition densities 
of 3no and 3.25no, respectively. 
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Figure 2: The density profiles of the star models considered in 
this work. The solid lines represent neutron star models with 
an APR equation of state, the dotted lines represent strange 
stars with a bag model equation of state and the dashed, dot- 
dot-dashed and dot-dashed lines represent hybrid star models 
with a large, medium and small quark matter core respec- 
tively. Thick lines represent 1.4 M0 stars and thin lines rep- 
resent massive 2 M© stars. In contrast to strange stars that 
are basicaUy homogeneous, stars that contain hadronic mat- 
ter have a very strong density dependence that extends over 
14 orders of magnitude reflected by the near zero segments in 
this plot. The very thin solid curve presents the maximum 
neutron star model for the APR equation of state ~ 2.2 Mq 
where hadronic direct Urea processes are allowed to the left 
of the dot. The dotted horizontal line denotes the density 
n = no/4 chosen as the beginning of the crust whose contri- 
bution is not taken into account in the damping time integrals 
below. 



B. R-mode profile 

The analysis of oscillation modes of compact stars re- 
quires the solution of the corresponding hydrodynamics 
equations [1]. Whereas the solution of the general rel- 
ativistic equations for static star models is straightfor- 
ward, the corresponding dynamic equations are more in- 
volved and explicit analytic expressions for the r-mode 
oscillation are only available in non-relativistic approxi- 
mation. Yet numerical analyses of the general relativis- 
tic equations show that except for very compact stars 
the relativistic corrections are moderate [34-36 . In the 
non-relativistic case the current state of the art is the 
comprehensive analysis It might therefore seem more 
consistent to employ Newtonian equations for the static 
star models as well, but we prefer to perform the neces- 
sary approximation for the oscillation at least around the 
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correct equilibrium configuration. In particular since we 
study heavy quark and hybrid stars where the difiierent 
approximations could differ. In particular in the large 
amplitude regime, a consistent general relativistic analy- 
sis would clearly be desirable. Furthermore, so far these 
analyses assume that the oscillation modes are solutions 
of an ideal fluid. For these modes the damping is then 
computed in a second step. 

R-modes are normal oscillations of rotating stars and 
correspondingly they require as a first step the solution 
of a uniformly rotating stellar model. Since neutron stars 
are cold, dense systems we assume a barotropic fluid 
where the pressure p is only a function of the energy 
density p. The hydrodynamic Euler equation for this 
spherical system, determining the enthalpy h and the 
equation for the gravitational potential $ of the star, 
have to be solved with appropriate boundary conditions 
[3]. To simplify the demanding analysis a slow rotation 
expansion is performed and, since the density fluctuation 
of the r-mode vanishes to leading order in the expansion, 
the computation of bulk viscosity damping times requires 
an expansion of /i, <1> and of the energy density p to next- 
to-leading order 



where the parameter k can likewise be expanded in Q 



X{r,cose) = Xo(r) +X2(r,cos( 



nGpo J 



where X stands for either of the three quantities, is the 
angular velocity of the rotation and po the average energy 
density of the corresponding non-rotating star. The main 
effect of the rotational corrections is a flattening of the 
star due to centrifugal forces. 

The next step is the search for eigenmodes of the rotat- 
ing star in a linear low amplitude approximation, e.g. for 
the conserved baryon number An <C n. They can com- 
pletely be described by the change in the gravitational 
potential and the hydrodynamical perturbation 
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These are likewise expanded to next-to-leading order in 
f2 in the form 



SX{r, cos 9) ^R^n^ [ dXo{r) + dX2{r,cose) 
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where SX stands again either for 5$ or SU and R is the 
radius of the static star. The potentials obey complicated 
differential equations with corresponding boundary con- 
ditions given in [3] . The boundary conditions require that 
the oscillation frequencies tJr in the rotating and Ui in the 
inertial frame are connected to the rotation frequency 
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We study classical r-modes which are a one parameter 
class of eigenmode solutions, I — m, that are to leading 
order determined by kq = 2/ (m + 1) and given by 



7r(m+l)'' (2to+1)! 
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in terms of associated Legendre polynomials P^^i and 
with the definition of the dimensionless amplitude a in- 
troduced in [6J. For the lowest r-mode that couples to 
gravitational waves with m = 2 the oscillation frequency 
in the inertial frame is given by uJi — —4/3 fi, correspond- 
ing to a counter-rotating flow. The leading order gravita- 
tional potential obeys a differential equation that is given 
in appendix [X| where its analytic solution is given in the 
special case that the star is of uniform density. In gen- 
eral it requires a numeric solution and then completely 
determines the r-mode to this order. 

At next-to-leading order there are two qualitatively dif- 
ferent effects: First the connection between oscillation 
and rotation frequency as described by the parameter k 
becomes non-linear via a non-vanishing value of K2 that 
has to be obtained numerically from a corresponding dif- 
ferential equation [3 . In flg. [3] the solution is shown for 
the star models discussed here. Whereas the corrections 
are small for quark stars, at large frequency they can be- 
come sizable for hadronic and hybrid stars. The second 
effect is the change of the potentials arising as solutions of 
rather involved partial differential equations. At next-to- 
leading order they feature nontrivial radial and angular 
dependences that are not described by simple power laws 
or spherical harmonics anymore. 

We are interested in the ampliflcation of these modes 
due to gravitational waves and their viscous damping 
which are described by the energy dissipation 

— = - U [UJ - mil) \dJm,n\ 



(fx {27^6(7"'' 5(7ab + CS(JS(7* 



(3) 



where rj and C are the shear and bulk viscosity, respec- 
tively. The fluctuations SJmm and Saab couple to gravita- 
tional waves and shear viscosity, respectively. The fluc- 
tuation (5(7 = V ■ 6v of an r-mode oscillation which is 
subject to dissipation via bulk viscosity reads 
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UJr = UJ — , UJi — UJr — mQ (2) 



where A denotes the inverse speed of sound 
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Figure 3: Connection of the oscillation frequency u)i of the 
r-mode in the inertial frame to the rotation frequency of the 
considered star models to next to leading order in the f2- 
expansion. The horizontal line shows the leading order result 
and the conventions for the other curves are the same as in 
fig- [21 
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evaluated at equilibrium and the dots represent several 
further terms of next-to-leading order in f2 given in [3] . 
The terms in the parenthesis depend on the expansion 
coefficients of the potentials and on the connection pa- 
rameter At, but are by virtue of the expansion independent 
of frequency. The fluctuation 5a is finally connected to 
the fiuctuations of energy Ap and baryon number density 
An via their continuity equations 
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Following in our numerical analysis we will consider 
the change of k at second order, but because of the in- 
volved numerics [3] , we will not take into account all ex- 
plicit second order terms in eq. Q. A standard approx- 
imation [B] is to replace the Lagrangian density fluctua- 
tion by the Eulerian one which corresponds to neglecting 
the terms denoted by the parenthesis in eq. Q. E.g. 
the density fluctuation for a m = 2 r-mode reads then to 
leading order 
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dp 
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in terms of the spherical harmonics Fg^. However, we 
will give general semi-analytic expressions below that are 
valid to full next to leading order. These show that the 
influence of the neglected terms on important aspects of 
the instability regions is rather mild. 



III. VISCOSITIES OF DENSE MATTER 
A. General expressions 

The bulk viscosity describes the local dissipation of 
energy in a fluid element which arises from the global 
oscillation mode of the star within a compression and 
rarefaction cycle. The integration over the whole star 
yields then the corresponding energy dissipation of the 
mode as will be discussed in the next section. Recently 
the bulk viscosity of large amplitude oscillations [I8_ has 
been studied in detail [ITj. There it was shown that 
large amplitude oscillations are in general considerably 
more strongly damped and this mechanism can thereby 
saturate unstable r-modes at flnite amplitudes as will be 
discuss in a companion article. In this article, however, 
we restrict ourselves to the subthermal regime /iA ^ T 
which determines whether small amplitude r-modes are 
initially unstable. 

The bulk viscosity is maximal when the external os- 
cillation frequency matches the time scale of the mi- 
croscopic interactions that establish equilibrium. The 
relevant interactions in the case of star oscillations are 
slow weak processes. The parametric form of the beta- 
equilibration rate is given by 




(7) 



where /za = X]j Mi ~ Pf quantity that is driven 

out of equilibrium due to the oscillations and its re- 
equilibration leads to the bulk viscosity. In the subther- 
mal regime, /ia ^ T, the non-linear terms in the rate of 
the equilibration processes can be neglected, leading to 
the general analytic result for the subthermal bulk vis- 
cosity [Tf] 
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2 Smaa; ^ J2 



(8) 



in terms of the reduced weak rate T and the strong sus- 
ceptibilities B and C 
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as well as the oscillation frequency w (in the rotating 
frame) and temperature T. It has a characteristic reso- 
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nant form and as long as the combination of susceptibil- 
ities / B does not vary too quickly with temperature, 
the sub-thermal viscosity has a maximum 
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max 
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and simple asymptotic limits 
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where / = BTT^ / uj can be identified as the feedback 
term in the differential equation that determines the 
damped oscillation |I2j. 

The shear viscosity arises from strong or electromag- 
netic interactions. In contrast to the bulk viscosity, the 
shear viscosity of dense matter is independent of the fre- 
quency of an external oscillation and approximately de- 
pends on temperature via a simple power law. Shear vis- 
cosity becomes large at low temperatures and therefore 
it is the dominant process for damping of the r- modes of 
cooler stars. Thereby, to leading order it can be param- 
eterized as 



(11) 



by simply factoring out the temperature dependence with 
exponent a. In general several processes can contribute 
so that the full shear viscosity can approximately be writ- 
ten as a sum of such power laws for the individual pro- 



B. Viscosities for the considered forms of matter 

In the following we discuss the viscosities for the phases 
of dense matter presented in the considered classes of 
compact stars. We note however, that the general results 
given here likewise apply to more complicated forms of 
matter like hyperonic and/or superfluid nuclear matter 
as well as various forms of superconducting quark matter. 



1. Bulk viscosity 

The general expression for the bulk viscosity eq. ([s]) 
depends on the weak rates. In the case of hadronic matter 
in the absence of hyperons the weak equilibration occurs 
via the Urea channel 



are necessary to satisfy energy momentum conservation. 
The latter represents a particular strong interaction ver- 
tex correction to the above process. However, from the 
point of view of the weak interaction these different pro- 
cesses belong to the same channel. Further, there is in 
principle also a corresponding process with muon instead 
of electrons. Since the phase space and thereby the rate 
of this process is strongly suppressed compared to the 
electron version at moderate densities we do not consider 
it here. 

In strange quark matter the dominant channel for beta 
equilibration is the non-leptonic flavor changing process 



s + u d + u 

whereas the corresponding quark Urea processes are 
parametrically suppressed in T/fi <JC 1. The weak pa- 
rameters r, S and X in the parameterization of the weak 
rate eq. ^ are given on the left panel of table II] 

Further, the inverse speed of sound A, eq. (5), as well 
as the strong susceptibilities B and C, eq. ([9 1, that pa- 
rameterize the deviation from chemical equilibrium are 
required. In order to make it easy to apply our general 
results to different forms of interacting hadronic matter, 
we implement the APR equation of state using the sim- 
ple parameterization employed in [41J to approximate the 
dependence of the energy per particle on the proton frac- 
tion x by a quadratic form 

E{n, x) = E,{n) + S{n) (1 - 2xf 

where and S are the corresponding energy for sym- 
metric matter and the symmetry energy for which we 
employ simple quadratic fits to the APR data. For com- 
parison we also consider the susceptibilities in the ap- 
proximation of a hadron gas which had been used in pre- 
vious analyses of the bulk viscosity [71 [33 HO] ■ Likewise 
we compute these quantities for the generic, phenomeno- 
logical form of the equation of state of interacting quark 
matter eq. Q. The resulting strong interaction param- 
eters describing the response of the different models are 
given in table ITTT] 



2. Shear viscosity 

In previous r-mode analyses the shear viscosity in 
hadronic matter has been approximated by the contri- 
bution from strong hadron-hadron-scattering using the 
fit in [42J to the standard low density (< no) data given 
in 113] 



n— >p + e + Ve , p + e n + Vf, 

There are two qualitatively different cases depending on 
whether the direct process |371 [35] is possible or only 
the modified version [38fl4D] where by-stander nucleons 



7?„ = 347pJT-2^ (12) 
cms 

where T is in units of Kelvin and p is given in g/cvo? . Ex- 
trapolating this fit to high densities relevant for neutron 
stars overestimates the viscosity. The new evaluation in 
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Weak process 


f [MeV*^"*'] 


5 


Xi 


X2 


X3 


quark non-leptonic 


6-59x10-1^(33^) 


2 


47r^ 








hadronic direct Urea 


5.24x10-'^ (1^)' 


4 


10 


1 

177r<l 





hadronic modified Urea 


4.68x10-^" 


6 


189 


21 


3 






3677r'* 


18357r6 



Strong/EM process 


77 [MeV(^+"'] 


a 


quark scattering 


b 14 


5 
3 


leptonic scattering 


1.40X10^^ /2L2i\T 

^""^ , — 


5 
3 


nn-scattering 


5.46 X 10'-' f " V 


2 



Table II: Left ■panel: Parameters of the general parameterization of the weak rate eq. ([7| for diflferent processes of particular 
forms of matter which determine the damping due to bulk viscosity. The coefficients Xi parameterize the non-linear dependence 
on the chemical potential fluctuation /xa arising in the suprathermal regime of the viscosity which is relevant for large amplitude 
r-modes studied in [16| but which are not relevant in this work. Right panel: Parameters arising in the parameterization eq. 
( |11[ | of the shear viscosity for different strong and electromagnetic interaction processes. The leptonic and quark scattering 
arises from a non-Fermi liquid enhancement due to unscreened magnetic interactions. 



hadronic matter 
hadronic gas 



quark matter (gas: c = 0) 



B 



3m^ 
(3ir2n) 3 



. ' (4(l-2x)Sy 
3(37r2) 3 n 3 



3 + 



27r^ 



(1^)H2 3(l-c)Ai2 



1 + 



12{1-c)m2 



c 



4(l~2.)(n^-f) 

(37r^n) 3 
6m 



3(l-c)Mq 



Table III: Strong interaction parameters, defined in eqs. Q and Q, describing the response of the particular form of matter. 
In the case of interacting hadronic matter a quadratic ansatz in the proton fraction x parameterized by the symmetry energy 
S is employed. The expressions for a hadron and quark gas are given to leading order in n/m% respectively next to leading 
order in nis/jJ.. 



|44| shows instead that due to a non-Fermi liquid en- 
hancement arising from the exchange of Landau-damped 
transverse photons, the main contribution to the shear 
viscosity of hadronic matter, at temperatures relevant to 
the spin-down evolution of the compact stars, comes from 
electron scattering and is given by 



4.26 X lO-'"'{xn)T^T-3. 



(13) 



where T is in Kelvin and the baryon number density n is 
in units of cm""^. In the calculation of rjf, we have only 
considered electron-electron and electron-proton scatter- 
ing and neglected the small effect of the muons to the 
shear viscosity. For densities larger than nuclear matter 
saturation density this electron contribution dominates 
over the hadronic one down to temperatures below 10*^ K. 
This region contains the part of the instability region that 
is relevant for the spin-down evolution of stars. There- 
fore in our main analysis we will completely neglect the 
hadronic component of the shear viscosity. However, we 
will compare with the previous form of the shear viscos- 
ity eq. ( 12 ) in order to discuss the effect of our improved 
analysis on the instability region^. 

In the case of ungapped quark matter, the shear vis- 
cosity is dominated by quark-quark scattering, and in the 



^ The leftmost part of the instability region below ~ 10*^ K only 
features an instability at very large frequencies. This part would 
only be relevant for old stars that are spun up by accretion, yet 
frequencies very close to the Kepler frequency are not reached 
via this mechanism anyway, cf. fig. ^ due to turbulent effects. 



limit of T <^ qjj, where is the Debye wave number, it 
is given by |J5j 



1 



Vq = 



407ra 



2N, 



(14) 



1.81, 



where = |^ is the QCD coupling constant, a 
Nq = 3 and /ig and T are in units of MeV. The temper- 
ature dependence arises again from a non-Fermi liquid 
enhancement of the quark interaction. These expressions 
yield the parameters in the parameterization of the shear 
viscosity eq. (11 1 as given on the right panel of table 



IV. R-MODE TIME SCALES 

A. General expressions 

The amplitude of the r-mode oscillations evolves with 
time dependence e'"*~*/'^, where uj is the real part of the 
frequency of the r-mode and 1 /t is the imaginary part of 
the frequency. The latter describes both the exponential 
rise of the r-mode driven by the Friedman-Schutz mecha- 
nism pS] and its decay due to viscous damping. We can 
decompose 1/r as 

1111 



T(f),T) ram TBin,T) ts{T) 

where tq, tb and are gravitational radiation, bulk vis- 
cosity and shear viscosity time scales, respectively. The 
lowest of these individual time scales determines if the 
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r-mode is unstable or damped. The damping time for 
the individual mechanisms is in general given by 



1 

n 



1 

'2E 



dE 
~dt 



(15) 



and requires both the total energy of the r-mode 



E 



1 



R 



2m+2 



dr 



(16) 



and the dissipated energy which is given by the corre- 
sponding part of eq. ([3| [3 . For instance the dissipated 
energy eq. ([3| due to the bulk viscosity reads 



dE\ 
'dF , 



K^n^ I d^x 



Ap 

p 



Ap 



where the dependence of the bulk viscosity on the con- 
served number density fluctuation amplitude has been 
expressed in terms of the conserved energy density. If 
the star consists of several shells consisting of different 
forms of matter with different viscosity, as is the case for 
hybrid stars, neutron stars with a high density core where 
direct Urea reactions are allowed, etc., the integral con- 
sists of partial integrals over the individual shells s and 
the inverse viscous damping times can be written as 



Therefore, the contribution of the different shells enters 
additively. We will in the following expressions sup- 
press the explicit label (s) but implicitly assume that 
the damping times consist of several contributions when 
different layers are present. With these expressions the 
individual r-mode time scales can be obtained. The time 
scale of the r-mode growth due to gravitational wave 
emission is given by [B] 



327r (m — 1)^™ / to 



1 _ 

TG ^ ((2to + 1)!!)^ V™+ 1 
with the radial integral constant 



2m-|-2 



(17) 



Jnn, 



1 



R 



p{r) r 



2m+2 



dr 



(18) 



1. Shear viscosity damping time 

The damping time of an r-mode with angular quantum 
number to due to shear viscosity is given by 



1 _ (to-1)(2to + 1) 



rir'^'^dr 







Using the parameterization eq. (Ill this can be written 
as 



(to - 1) (2to + 1) S„,K 



1 

in terms of the dimensionless constant 

n-Ro 



3 + <T 

QCD 



R 



Sm = 



1 



^^QCD -JRi 



fir^"'dr 



(19) 



(20) 



where Ri and Ro are the inner and outer radius of the 
corresponding shell, if there are several ones. This con- 
stant contains the complete dependence on the particular 
microscopic processes. To make this quantity dimension- 
less the generic scale Aqcd has been introduced that is 
chosen as Aqcd = 1 GeV for the numeric values. For the 
fundamental m — 2 r-mode the parameter S is given for 
the different star models considered here in table HVl 



2. Bulk viscosity damping time 
The bulk viscosity damping time is given by 



1 



tb a'^JmMR^ 





An 




An 




1 d^x 




•<( 








n 




n 





(21) 



where the bulk viscosity is in general a function of the 
amphtude |T71|TH|. As has been noted before the strongly 
enhanced damping can provide a mechanism for the satu- 
ration of the r-mode as will be discussed elsewhere. Here 
we restrict ourselves to a study of the initial instability of 
small amplitude r-modes and in the subthermal regime 
the viscosity eq. ^ is independent of the r-mode am- 
plitude. In this case only the angular integral over the 
density fluctuation enters and it is useful to define the 
angular averaged form 



S^{r) 



m 



1 



m+l 



(to + 1)''^ (2to + 3) 
4to 



dn\6a\ 



l.o. 



(- 

\R 



<5$o (22) 
which reduces to the second line to leading order in the Cl- 



expansion, where the angular integral in eq. ( 21 1 is trivial 
due to the normalization of the spherical harmonics. The 
damping time in the subthermal regime is then given by 



16to 
(2to + 3) (to - 



1? K Jm M 



(23) 



10 



where the dependence on all local quantities, like the 
equation of state, the weak rate, the density dependence 
of the particular star model and its r-mode profile is con- 
tained in the function 



dr 



iSnr)f (24) 



depending on a single external parameter and which has 
to be determined numerically for a given star model. In 
the asymptotic limits the damping time simplifies to 



16to 



/«i (2to + 3) (m + 1)^ k2 



16m 



/»i (2m + 3)(m + l)^ 
with the dimensionless constants 



(25) 



(26) 



a4 i-Ro 



w„ 



/?3/\9-<5 , 



ffSA^ A''^^ Id 



drr'A'C'T {5Y.{r)Y 



(27) 



Tip 



dr r' 



mr)f (28) 



Since the bulk viscosity originates from weak interactions 
here the second normalization scale Aew is used with 
a generic value Aew = 100 GeV. These normalization 
scales are only introduced to obtain dimensionless con- 
stants of order one and drop out of the final results for 
the damping times. 

Finally, we want to stress that in these expressions for 
the r-mode time scales the complete local dependence 
on the star profile, the r-mode oscillation and the micro- 
scopic damping processes is contained in the few constant 
J, S, V and W, but the dependence on the global pa- 
rameters of the r-mode evolution fl and T is entirely 
explicit. These constants include in particular also the 
complete dependence on the non-trivial radial and angu- 
lar dependence of the full next-to-leading order expres- 
sion for the r-mode |3j. Since to our knowledge results for 
the bulk viscosity in the crust of neutron stars are not 
yet available, we neglect the damping of the crust in our 
numerical analysis. It has, however, been argued that 
the shear viscosity of the crust could be crucial [TS] and 
a more detailed study of this issue is definitely desirable. 



B. Results for the considered star models 

Before we discuss the damping times for particular star 
models let us point out a few generic properties of the 
damping times that follow directly from the general ex- 
pressions eqs. (17 1, (19 1 and (23 1. The gravitational 



time scale is independent of temperature and decreases 



strongly with frequency. The shear viscosity damping 
time increases with temperature and is independent of 
frequency. Because of the resonant behavior of the bulk 
viscosity, the corresponding damping time decreases with 
temperature, has a minimum and increases again at large 
temperatures. Furthermore, it also decreases with fre- 
quency, but slower than the gravitational time scale. 

A numeric solution requires the integration over the 
star profiles to obtain the constants J, S and the func- 
tion T^, respectively the constants V and W describ- 
ing its asymptotic behavior. The constants are given for 
the fundamental m — 2 r-mode of the various star mod- 
els discussed in this work, and where applicable also for 
their different shells, in table \rv\ For the bulk viscosity 
of a 1.4 Mq neutron star we compute the APR result em- 
ploying the proper susceptibilities for interacting matter, 
as well the result when the susceptibilities are evaluated 
in the idealized case of a hadronic gas, as has been done 
previously |40ll42j . As can be seen the parameters in the 
interacting case are more than an order of magnitude 
larger, owing to the larger viscosity [17^, which leads to a 
correspondingly smaller damping time. In the following 
we will consider only the proper interacting form unless 
otherwise noted. 

In the left panel of fig. [4] the numeric solution for the 
different time scales is shown as a function of temperature 
for 1.4 M0 stars of the different classes considered here 
rotating at their Kepler frequency flu- Shown are grav- 
itational time scales (horizontal lines), the shear viscos- 
ity damping times (monotonically increasing curves) and 
bulk viscosity damping times (non-monotonic curves). 
The solid lines denote the neutron star model, the dotted 
lines a strange star and the dashed and dot-dashed lines 
show hybrid star models with a small and large quark 
core, respectively. Whereas damping due to shear viscos- 
ity dominates for strange and hybrid stars for T < 10^ 
K and for neutron stars even for T < 10^ K, the bulk 
viscosity damping time of the neutron and the strange 
star feature the generic resonant form, where the min- 
ima are around 10^ K and 10^^ K, respectively, and the 
higher gradients in case of the neutron star arise from 
the higher power of 6 (tab. |ll|. It is clear that as a func- 
tion of the core size the family of hybrid curves interpo- 
lates continuously between the two uniform star models. 
Correspondingly, for hybrid stars the contribution of the 
quark core dominates at low temperature, whereas the 
hadronic shell dominates at higher temperatures, which 
leads to a curve with two minima. For all star models at 
ilfcthe gravitational time scale is the lowest over a range 
of temperatures. 

The right panel of fig. |4] shows the same plot at the 
lower frequency fix / 4. As is clear from the general dis- 
cussion the shear curves are unchanged whereas the grav- 
itational curves and the bulk curves are shifted upwards 
compared to the left panel of fig. [4] Since the gravi- 
tational time scale moves upwards faster than the bulk 
viscosity curve, the instability region shrinks from both 
sides as the frequency is lowered. For hybrid stars with 



11 



^ 10' 
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T[K] 



Figure 4: Damping times of the different 1.4 star models discussed in this work. Shown are a hadronic star with modified 
Urea processes (solid), hybrid stars with a small (dot-dashed) and large (dashed) quark core as well as a strange star (dotted). 
The horizontal curves give the time scale re associated to the growth of the mode due to gravitational wave emission. The 
monotonically increasing curves show the damping time ts due to shear viscosity and the non-monotonic curves the damping 
time tb due to bulk viscosity. Left panel: Stars rotating at their Kepler frequency fix- Right panel: Same for stars rotating at 



star model 



shell 



J 



S 



V 



w 



NS 1.4 M( 







1.81 X 10" 



NS 1.4 Mq gas 



7.68 X 10'^ 
4.32 X 10"*^ 



NS 2.0 M( 







2.05 X 10" 



2.25 X 10" 



1.31 X 10"^ 
1.28 X 10""^ 



1.61 X lO"*^ 



3.49 X 10' 



371 



5.05 X 10""* 



1.16 X 10"^ 



1.76 X 10"'^ 



3.70 X lO'' 



226 



1.16 X lO"-** 
9.34 X 10"* 



1.72 X 10" 



4.19 X 10' 



368 



NS 2.21 Mp 



d.U. core 
m.U. core 



2.02 X 10" 



7.11 X 10""^ 
1.55 X 10"® 



2.03 X lO'' 



493 



SS eq. 



SS 1.4 M, 



all 



287r 



© 



3.08 X 10" 



, 14/3 

3.49 X lO"*^ 



47rlAt 



3.53 X 10" 



0.191 



eq. 



34 



eq. 



33 



eq. 



36 



eq. 



35 



7.86 X 10'^ 



1020 



1.01 X 10" 



8340 



SS 2.0 M, 



© 



2.65 X 10" 



4.45 X 10"*^ 
2.19 X 10"*^ 



3.38 X 10"^° 



0.157 



..58 X 10'^ 



955 



;.07 X 10* 



6600 



HS 1.4 M0 L quark core 



1.93 X 10" 



1.38 X 10"^" 



1.03 X 10" 



9.87 X 10' 



1170 



5.83 X 10** 



7270* 



hadr. core 



6.72 X 10"' 
2.88 X lO"'' 
3.29 X lO"*^ 



1.40 X 10" 



1.51 X lO"'^ 



X 10' 



1040 



6.85 X 10~^ 
1.56 X 10"® 



HS 1.4 Mq S 



quark core 
hadr. core 



X 10" 



7.96 X 10"" 
1.25 X 10"^ 



1.01 X 10^ 
3.61 X 10' 



1070 

394'' 



4.35 X lO*' 



3690* 



5.25 X 10"® 



3.76 X 10"" 
1.07 X 10"^ 



2.34 X 10"^ 
1.32 X 10"® 



HS 2.0 Mq quark core 



2.00 X 10" 



7.73 X 10' 



1080 



4.77 X 10* 



6310* 



hadr. core 



5.24 X 10"® 



7.65 X lO'' 



925f 



Table IV: Radial integral parameters ajid_ characteristic points of the i nsta bilit y r egio n of am = 2 r-mode for the star models 
considered in this work. The constant J, S, V and W are given by eqs. ( 18 1, ( 20 1, ( 27 1 and ( 28 1 using the generic normalization 



scales Aqcd = 1 GeV and Aew = 100 GeV. The temperatures and frequencies are obtained with the analytic expressions for 
the minima eqs. (B2|, (B3l, (B81 and (B9l and for the maxima eqs. (B5l and (B61. The expressions for a generic strange 
star (or quark core) in terms of the parameters of the quark model equation of state eq. |l]), using the constants f and fj 
defined in eq. (32 1 are given as well. (*) These values deviate significantly from the actual results due to the inappropriate 
approximation of constant radial profiles, whereas this idealization entirely fails for hadronic parts. (|) The second minimum 
arises from the competition of the bulk viscosity damping in the quark and the hadronic shell, (o) The second minimum arises 
from the competition of bulk and shear viscosity damping in the hadronic shell. 
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a large quark core the gravitational time scale can move 
above the minimum of the bulk viscosity arising from the 
quark core as the frequency is lowered and the instability 
region is divided by a stability window. Yet, for a hybrid 
star with a sufficiently small quark core the minimum 
can move above the shear curve before the gravitational 
time scale can overtake it. In this case there would be no 
signature of the quark core and the instability regions of 
such a hybrid star and of the corresponding neutron star 
would look basically indistinguishable'*. This is nearly 
the case for the small core hybrid star model denoted by 
the dot-dashed curve as will be discussed below. 



numeric solution of eq. (29) is straightforward, these an- 



alytic expressions reveal the general dependence on the 
various unknown model parameters entering the r-mode 
analysis and therefore provide a measure for the uncer- 
tainty of these results in a situation where most proper- 
ties of dense strongly interacting are basically unknown. 
Here we give expressions for the fundamental to = 2 r- 
mode and the classes of stars discussed in this work and 
defer the general results to appendix [B| 



1. Low temperature boundary 



V. R-MODE INSTABILITY REGIONS 

The boundary of the instability region is given by the 
condition \tq\ = \tv\, where Ty is the viscous damping 
time. This reads in terms of the individual contributions 



As is clear from the analytic expressions for the damp- 
ing times and confirmed by fig. [4] at low temperatures 
shear viscosity damping dominates. Furthermore, in the 
case of hybrid stars the damping due to hadronic shear 
viscosity dominates over the quark core. Comparing it 
with the gravitational time scale yields 




= 



(29) 



Since at the boundary the time scale of the gravitational 
instability is identical to the viscous time scale, the r- 
mode amplitude neither increases nor decreases. If the 
viscous damping time is shorter modes induced by exter- 
nal perturbations will quickly be damped away, whereas 
in the opposite case they are unstable and will initially 
grow. As will be shown in a second article the increase in 
the bulk viscosity at large amplitude can eventually sat- 
urate the r-mode, but here we will limit ourselves to the 
small amplitude regime and analyze the regions where 
small amplitude modes are unstable. In general eq. (29) 



has no analytic solution and has to be solved numeri- 
cally. However there are various limiting cases for which 
analytic solutions exist and we will study them below. 



A. Analytic expressions 



The boundary of the instability region eq. (29) can 



generally not be found analytically due to the occurrence 
of several terms with non-trivial temperature and fre- 
quency dependencies. Since the viscous damping times 
vary extremely strongly with temperature and frequency 
(fig. |4]) in general there is over nearly the whole parameter 
space one component that clearly dominates the others in 
the sum. In such a case the equation can easily be solved 
analytically and yields analytic expressions for the dif- 
ferent segments of the instability region. Although the 



* The fact that the volume of the neutron shell is slightly smaller is 
not significant for the damping times that vary over many orders 
of magnitude, so that the dashed curves in fig. |4] are invisible 
underneath the solid neutron star curves. 



-> 1.12 , , , T- 



(NS andSS) 



where the irrational power arises from the Landau damp- 
ing of the corresponding interactions that induce shear 
viscosity. Interestingly this expression depends only 
mildly on the constant S that encodes the microscopic 
interaction. This segment of the boundary is important 
for old stars in binary systems that are spun up by ac- 
cretion from a companion star and enter the instability 
region from below 



2. Minimum of the instability region 

Since the shear viscosity damping time monotonically 
increases with temperature whereas the bulk viscosity 
damping time monotonically decreases at low tempera- 
ture there is a minimum of the instability boundary given 
by tb = ts- Due to the non-linear frequency dependence 
of the second order damping times an exact analytic so- 
lution is not possible. However, since these minima are 
located at frequencies far below the Kepler frequency it 
is according to fig. [3] a very good approximation to use 
the leading order frequency connection k = kq in which 
case an analytic solution is possible. For neutron stars it 
reads 



{NS) 



{NS) 



1.06- 



1.89 



49 ""23 23 23 ^ 

i? 128 J 64 G 128 M 64 A|\;j^ 

K 64 _K64 



(30) 



(31) 



Note, the appearance of surprisingly low powers in these 
expressions. In particular, due to the arising 5/128 power 
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a change of V by an order of magnitude results only in 
a mild deviation of the minimum frequency of less than 
10%, and even a very drastic change by three orders of 
magnitude does not change the result by more than 30%. 
Whereas the viscosity constants eqs. (27 1 and (28 1 can 



vary by many orders of magnitude for different classes 
of stars with different transport processes as can be seen 
from tab. |IV[ they are generically of similar order of 
magnitude within a given class due to the identical para- 
metric dependence on the microscopic physics. Recall, 
that the minimum of the instability region is of particu- 
lar importance for the r-mode analysis since it determines 
to what frequency r-modes can spin down a star. It is 
needless to say that such an insensitivity of the minimum 
frequency on the microscopic physics is more than wel- 
come in the present situation where there are still huge 
uncertainties on the underlying equation of state and the 
transport coefficients of dense matter. This presents one 
of the main results of this article. 

In the general case studied in appendix |B] the 
{2m {6 + a) + 2S) /a-th root of the bulk viscosity con- 
stants Vm arises in the expression for the minimum fre- 
quency eq. ( B2 I . For short-range (Fermi liquid) inter- 



actions the shear viscosity exponent is generically a — 2 
whereas it reduces to cr — 5/3 when long-ranged, only 
Landau-damped (non-Fermi liquid) gauge interactions 
are present. Since the rate of the weak interactions van- 
ishes in equilibrium and requires phase space both for ini- 
tial and final state particles in general S > 2. A similar 
argument should hold for processes mediated by Gold- 
stone bosons in color superconducting phases. Therefore 
the above root is for all multipoles and corresponding 
processes higher than 11 so that change of the viscosity 
constants Vm by an order of magnitude still changes the 
minimum only by at most ~ 20%, making it very in- 
sensitive to them. This extends the finding obtained in 
an explicit comparison of various different neutron star 
models in [B] to general forms of matter and damping pro- 
cesses. This insensitivity is particularly interesting since 
these constants contain the complete dependence on the 
detailed second order r-mode profiles, so that the min- 
imum frequency is hardly affected by the second order 
effects. This observation had already been made for the 
particular neutron star model studied in [3]. A similar 
statement holds for the dependence on the shear viscos- 
ity constants Sm , where the corresponding root is at least 
of 6th order reached in the limit 6 ^ a. The tempera- 
ture of the minimum is more sensitive to these constants, 
as likewise observed in [31 |B] , since it only involves the 
(cr + (5)-th root. 

Finally we give the explicit expression for the minimum 
frequency of strange stars with the general quark model 
equation of state eq. ([T|. Here the parameter dependence 
of the quark matter viscosity coefficients given in table 
ITT] can be factored out according to 



where F and fj are pure constants, so that the minimum 
is located at 



(SS) 



t: 



{SS) 



2.23- 



2.78 



-^5 3 

G5B(1 -c)AaFi?iM55 
fjTiG^^f (1 - c)^ Mn 



(33) 
(34) 



Assuming that eq. ([T]) gives an estimate for the uncer- 
tainty in the unknown quark matter equation of state and 
estimating the uncertainty in the quadratic and quartic 
parameters mg and 1 — c each by a factor of two and 
that in the strong coupling as generously by an order of 
magnitude, the minimum could vary here by roughly a 
factor of two owing to the uncertainty in the unknown 
microscopic dynamics. As noted before the quark model 
equation of state is also valid for color superconducting 
phases [32J. In this case the gap reduces the parame- 
ter nig which would lower the minimum frequency and 
enlarge the instability region. 



3. Intermediate boundary 

Above the minimum the damping is dominated by the 
bulk viscosity and below its maximum it can be approxi- 
mated by the asymptotic low temperature form eq. ( 25 I . 



The semi-analytic result to next-to leading order in is 
given in the appendix, but because of the small effect 
of the next-to-leading order corrections of the oscillation 
frequency on the instability regions, observed above, we 
can simply neglect these and obtain 
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This part of the instability boundary is important for 
the spin down of young compact stars since it determines 
where they hit the instability region during their initial 
fast cooling phase. 



4- Maximum of the instability region for strange stars 

Due to the resonant form of the bulk viscosity eq. ^ 
whose maximum eq. ( 10 ) translates into a minimum 



(32) 



of the corresponding damping time, the instability re- 
gion features a maximum if the corresponding frequency 
is below the Kepler frequency or otherwise splits into 
two parts. Since there is no analytic expression for the 
damping time in the vicinity of the maximum there is no 
general expression for the maximum. Yet, for hadronic 
matter the corresponding maximum frequency is usually 
above the Kepler frequency anyway so it is not relevant 
for physical applications. In contrast, for strange stars. 
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where the maximum is generally below the Kepler fre- 
quency and important for the r-mode analysis, the den- 
sity profile varies only mildly and can be approximated 
by a uniform density. In this case a general analytic 
expression for the maximum of the instability region is 
given in eqs. (B5 B6l in appendix p] and reduces using 



the susceptibilities in tab. |III|to 



riL^fj « 0.434- 



mi i?3 



t£^) « 0.210 



(l-c)5GsAfi 
(1 - c)5 mfiiB 



(35) 
(36) 



where the tiny second order corrections to the oscillation 
frequency, cf. fig. [3j were neglected in this case. Note 
that analogous to the expression for the maximum of 
the viscosity of quark matter given in [T7] the chemical 
potential drops out and therefore there is no ambiguity 
where to evaluate the susceptibilities in case the density 
distribution of the considered star model is not entirely 
constant. This gives for the dimensionless ratio 



n 



iSS) 



1.03- 



^3 p 6 



K 



(37) 



where miso. Mi. 4 and Rio are the effective strange quark 
mass in units of 150 MeV, the stars mass in units of 
IAMq and the radius in units of 10 km, respectively. 

Whereas the dependence on the parameter c, which 
in the perturbative regime is positive and thereby in- 
creases the maximum frequency, is comparatively mild, 
the dependence on the effective strange quark mass is sig- 
nificant and can within the probable uncertainty region 
100 MeV < nig < 200 MeV strongly change the position 
of the maximum frequency from values considerably be- 
low to values above the Kepler frequency. For supercon- 
ducting matter the gap can reduce the effective param- 
eter rus- Whereas the maximum frequency of massive 
stars can be significantly lower, the radii of sufficiently 
massive strange stars do not vary much with mass along 
a mass-radius curve so that the dependence on the ra- 
dius is mild despite the arising high power. To judge 
the implicit dependence of the masses and radii on the 
equation of state it is useful to recall that the solution 
of the TOV equations for a uniform density star exhibits 
scaling with the bag constant R/Rq,M/Mo ^ ^Bq/B 
[47| . Despite the uniform density approximation, the an- 
alytic result eq 



(37 



for the maximum, given in tab. IV 



agrees nicely with the numeric results for the considered 
strange star models. The above expression can also be 
applied to the corresponding analysis of Jaikumar, et. al. 
[7]. Using their lower strange quark mass mg = 100 MeV 
as well as their slightly smaller radius and taking into 
account that these authors erroneously used the angular 
velocity in the inertial frame to evaluate the bulk vis- 
cosity yields Vlmaxl^K ~ 0.47 in very good agreement 



with the plot of the exclusion region shown in their fig. 4. 
When employing the proper angular velocity in the ro- 
tating frame, eq. (37 1, gives instead the corrected result 
O^'^'^) /o 



0.59. 



The analytic results in appendix |B] to some extent also 
apply to the maxima of the instability regions of hybrid 
stars but in this case the agreement is less precise since 
the density profile of their quark core features a more 
pronounced radial dependence. 



5. High temperature boundary 

For completeness we also give the high temperature 
part of the boundary. 



T2>T„ 




Due to the resonant behavior of the bulk viscosity which 
decreases again at large temperatures, the r-mode is un- 
stable for all forms of matter at large temperatures above 
the corresponding boundary of the shell with the largest 
resonance temperature eq. ( 10 1 . This part could be less 



relevant for the r-mode evolution since stars cool very fast 
initially and could leave this region before the r-mode can 
develop. 



6. General properties of the analytic analysis 

Finally let us note a few generic properties of these 
analytic results. Since the mass appears in the denom- 
inator in all these expressions, an increase of the mass 
of the star increases the instability region and moves it 
uniformly downward and slightly to the right in a T — fi- 
plot. The dependence on the radius is less uniform, but 
in general these expressions are less sensitive to the ra- 
dius than to the mass of the star. Furthermore for a 
given equation of state the radius of a star varies very 
little for masses of 1 to 2Mq. As already noted, the de- 
pendence of the minima of the instability region on the 
viscosity parameters is extremely mild. Due to the aris- 
ing small powers 0(1/10) even a large variation of the 
viscosity parameters by three orders of magnitude within 
the set of possible equations of state for a given class of 
stars, like e.g. pure neutron stars would not change the 
minimum by more than a factor of two. In contrast the 
dependence of the maxima on the viscosity parameters is 
far more pronounced and as discussed above changes by 
more than a factor of two are here easily possible. 



B. Numeric results 

Let us now discuss the numeric results for the exclusion 
regions and compare them to the semi-analytic expres- 
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sions obtained in the last section in order to assess the 
quality of the latter. Fig. [5] compares the expressions for 
particular segments of the instability region and for the 
extrema to the numeric solution. The left panel shows 
the results for a standard IAMq neutron star model 
(solid) and a corresponding strange star model (dotted) . 
The thick curves show the numeric results and the thin 
curves the analytic approximations and the dots denote 
the minima. Since we neglect the contribution from the 
neutron star crust there is only a single shell in each 
case. As can be seen, except for the regions around the 
extrema, the analytic expressions for the different seg- 
ments approximate the boundary of the instability re- 
gion extremely well and are mostly hidden underneath 
the numeric curves. Yet, the extrema are in turn well 
described by the corresponding analytic expressions. As 
noted in m [71 [TT] , due to the resonant form of the bulk 
viscosity there is a stability window around 10^ K where 
strange stars are not unstable against r-modes up to large 
frequencies. For the considered strange star model the 
maximum of the instability region is above the Kepler 
frequency so that two separate regions appear, but as 
discussed before the position of the maximum depends 
on the particular microscopic parameters. Due to the 
same qualitative resonant structure of the bulk viscosity 
there is also a second instability region at high tempera- 
tures > 10^^ K for neutron stars that has previously been 
neglected due to the employed low temperature approxi- 
mation to the hadronic bulk viscosity f40l |42] . Although 
the r-mode is initially unstable the naive expectation is 
that this should be irrelevant for the spin-down evolution 
since the star cools extremely fast at such high temper- 
atures and might leave this region before r-modes can 
develop. Yet, since the interior of the star is initially 
opaque to neutrinos and the cooling is delayed this is not 
entirely clear and requires further study. 

On the right panel of fig. |5] the comparison for the 
IAMq hybrid star model is shown. Here there are 
contributions from the different layers, but strikingly 
the analytic approximation works even in this more 
complicated case. This plot shows nicely the generic 
structure of the boundary of the instability region, dis- 
cussed in the previous subsection, when several shells 
are present that feature qualitatively different damping 
mechanisms. Around its resonant temperature the bulk 
viscosity damping mechanism of a given shell in general 
clearly dominates over those of the other shells. If the 
resonant temperatures of the different bulk viscosities are 
sufficiently separated there are several maxima, that de- 
fine corresponding stability windows. Below each of these 
maxima there is a minimum where a dominant mecha- 
nism is replaced by the next. However, in case two reso- 
nant temperatures are too close or a shell is too small to 
have a sizable impact, individual stability windows can 
be fully or partly washed out. We discuss such cases 
below. 

Fig. [6] shows the exclusion region for the standard 
IAMq neutron star model compared to approximations 



used previously in the literature. The solid curve shows 
our new result with shear viscosity due to the dominant 
Landau-damped lepton scattering |33] and bulk viscosity 
due to modified Urea processes based on susceptibilities 
for interacting matter [T7]. The dotted curve shows the 
exclusion region when using the result for the shear vis- 
cosity from the fit in [42' to the standard low density 
(< no) data from hadron-hadron-scattering given in |43] . 
Extrapolating this fit to high densities relevant for neu- 
tron stars overestimates the viscosity compared to its ac- 
tual, subleading size [33], leading to the smaller exclusion 
region. The dashed curve shows the result when employ- 
ing the previously used expression for the bulk viscosity 
|40[ |42] which employs susceptibilities in the approxima- 
tion of an ideal hadron gas, see table It is again the 
insensitivity to the viscosity parameter V discussed in 
the previous subsection that is responsible for the fact 
that these corrections, which change the damping time 
by more than an order of magnitude, have only such a 
mild effect on the exclusion region and the minimum eq. 
( 30 1 . It is interesting that the effect of the interactions is 



opposite to that of the 2nd order 57-corrections [3]. The 
combined curve formed by the dotted and the dashed 
segment is the exclusion region that had previously been 
studied in the literature. As can be seen the combined 
effect of both corrections is to move the instability re- 
gion to lower temperatures so that it extends to roughly 
10^ K. This is relevant for old stars at low temperatures 
that are spun up by accretion since the r-mode becomes 
unstable already at lower frequencies. 

The left panel of fig. [7] compares the instability regions 
of the different 1.4 Mq star models considered here and 
also includes the analytic results for the extrema. The 
minima in these plots corresponds to the maxima where 
two damping mechanisms cross in fig. |4] whereas the max- 
ima in fig. [7] correspond to the minima in fig. [4] As can 
be seen the instability regions of the hybrid stars interpo- 
late between the neutron star and the strange star curve 
as the size of the quark core increases, even though all 
these curves are based on rather different equations of 
state distinguished by the interaction parameter c, see 
table [ij The right panel of fig. [7] shows the instability re- 
gions for the 2 Mq models discussed in this work. As has 
already been observed as a generic feature of the ana- 
lytic analysis, the r-mode instability in heavy stars is en- 
hanced and the boundary of the instability region moves 
slightly to lower frequencies. As predicted by the mass 
dependences of the analytic expressions this is most pro- 
nounced in the vicinity of the maxima and milder in the 
vicinity of the minima. The approximate analytic results 
for the maxima deviate slightly from the maxima of the 
numerical curves. The reason for this is that the uniform 
density approximation is not justified in this case due to 
the strong radial dependence of the density profiles, see 
fig-i 

Fig. |8] shows the instability region for an ultra-heavy 
neutron star M ss 2.2 Mq denoted by the solid curve. 
Here, the instability region has grown to the point that 
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Figure 5: Comparison of the numeric results for the instabihty region (thick hnes) with the various approximate semi-analytic 
expressions (thin lines and blobs) presented in the text for the different classes of compacts stars. The dots present the results 
for the extrema and the thin lines which are valid away from the extrema represent the corresponding analytic results taking 
into account only the contribution from the dominant process and shell in the respective region. Left panel: Instability regions 
for the 1.4 neutron star (solid) and strange star model (dotted). Right panel: Same for the 1.4 A/0 hybrid star model. The 
thin dot-dashed curve on the right panel which deviates slightly close to the maximum is the leading order result without the 
frequency corrections in fig. [3] whereas the corresponding curves are indistinguishable on the left panel. Note that here in the 
following where the ratio Q,/Q,k is plotted this ratio is taken with the respective Kepler angular frequency for each star model, 
see table IT] so that the same value of Q/Qk corresponds to a different value of Q, for different curves. 



the lower and the upper part are about to merge. For this 
star the densities that are reached are high enough that 
direct Urea processes are kinematically allowed within 
an inner core with radius of roughly R/2. The enhanced 
damping due to direct Urea bulk viscosity in the inner 
core leads to a notch at the right side of the instabil- 
ity region, as also found in [15]. Yet, the modification is 
rather mild in this case: the dashed curve shows the same 
star model when direct Urea processes are artificially sup- 
pressed. In the opposite extreme when direct Urea pro- 
cesses are artificially allowed in the entire core, shown 
by the dotted curve, the enhanced damping does lead to 
a significant change of the instability region. Moreover, 
the frequency of the minimum of the instability region 
increases. Due to the resonant form of the bulk viscos- 
ity, however, the instability region does not uniformly 
shrink over the whole temperature range where bulk vis- 
cosity dominates, but the stability window moves to lower 
temperatures leaving the r-mode unstable at higher tem- 
peratures where it is otherwise stable in the presence of 
modified Urea processes. The temperature scales of the 
corresponding stability windows agree with the tempera- 
ture scales of the resonant maximum of the bulk viscosity, 
see [IT]. Direct Urea processes are very sensitive to the 
proton fraction of dense matter. Whereas the required 



fraction is roughly 14% in the case of the APR equation 
of state, reached at relatively high densities > 5no, this 
could be different for other equations of state. In a case 
where the direct Urea core is larger but a modified Urea 
shell is still present so that both of them have a sizeable 
volume fraction, their combined damping could lead to a 
larger stability window, yet still at parametrically larger 
temperatures than the stability window of strange stars. 

Fig. [9] shows the result for the higher multipoles for 
the different IAMq star models. As can be seen, the 
right boundaries of the instability regions for the higher 
multipoles are extremely close to the boundary of the 
fundamental m = 2 mode. Correspondingly these modes 
could easily be excited once the evolution of a star enters 
the instability region. Each mode that is triggered would 
lead to a further enhancement of the spindown of the 
star and could change the evolution. This could be par- 
ticularly relevant for young neutron stars that enter the 
instability region at high temperatures when the cooling 
is still fast but there are no strong reheating mechanisms 
which could prevent the star from substantially penetrat- 
ing the instability region. 

Fig. [To] finally compares our numeric results to astro- 
physical data. Although the rotation frequencies of pul- 
sars are known to high accuracy, the relevant core tem- 
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Figure 6: Modification of the instability region of a 1.4 
neutron star due to improved approximations for the micro- 
scopic transport properties. The solid curve shows the stan- 
dard neutron star model with shear viscosity due to the dom- 
inant Landau-damped lepton scattering and bulk viscosity 
based on the proper susceptibilities for interacting matter (as 
all other neutron star results in this work) . The dashed curve 
shows the result when using the fit given in |42| to the density 
data < no for the shear viscosity from hadron-hadron scatter- 
ing obtained in |43) . The dotted line shows the result when 
employing the previously used expression neglecting interac- 
tions to the susceptibilities contributing to the bulk viscosity 

EOllli]. 



peratures are not known for most pulsars and involve 
large uncertainties even for those stars where estimates 
are available. The reason is that the temperature of the 
surface is indirectly inferred from the data and it is neces- 
sary to abstract from it the core temperature via models. 
Such an analysis has been performed for two low mass X- 
ray binary systems Aql X-1 and SAX J1808.4-3658 |49j, 
where the pulsars accrete matter from companion stars. 
The horizontal lines give an optimistic estimate for the 
uncertainty of the temperature. As can be seen these 
data points are well within the instability region of neu- 
tron stars. This holds in particular for the faster one of 
the two pulsars, Aql X-1, which is right in the middle. 
Our analysis shows this statement cannot be undermined 
by the unknown equation of state. So if r-modes would 
spin down stars so fast that the stars would leave the 
unstable region on time scales short compared to astro- 
physical ones these stars could not be pure neutron stars 
and would require some form of enhanced damping. Yet, 
in the case of strange stars and hybrid stars with a suf- 
ficiently large quark core, the data points are close to 



the boundary of the corresponding stability window and 
would be compatible in such a scenario^. As our results 
show, hybrid stars with a small quark core or neutron 
stars with enhanced hadronic direct Urea damping, in 
contrast, would not be sufficient. 



VI. CONCLUSIONS 

Using general expressions for the viscosities of dense 
matter we have derived semi-analytic results for the 
damping times of small amplitude r-mode oscillations 
and the boundary of the instability region. Our results 
show that the boundary of the instability region and in 
particular its minimum, which determines to what ex- 
tent r-modes can spin down a fast star, are extremely 
insensitive to the quantitative details of the microscopic 
interactions that induce viscous damping. However, the 
instability regions can nevertheless effectively discrimi- 
nate between qualitatively different classes of stars. In 
particular strange stars and hybrid stars with sufficiently 
large quark cores feature a stability window that can- 
not be reproduced with standard neutron stars without 
some admixture of exotic matter that provides enhanced 
damping. We find that the presence of some form of 
exotic matter does not automatically lead to a stability 
window since the instability region of hybrid stars with 
a sufficiently small quark core is almost indistinguish- 
able from that of a neutron star. Similarly, the presence 
of neutron matter with direct Urea interactions will in 
most cases not considerably change the instability region. 
However, due to the demonstrated insensitivity of the in- 
stability regions to quantitative microscopic details, the 
clear determination of a very fast pulsar with ft > 300 
Hz in the temperature range around 10^ K could provide 
a convincing signature for some exotic form of matter. 
What remains to be shown in order to transform this 
into a strict signature is that the crust does not domi- 
nate the damping, and that r-modes do not saturate at 
amplitudes that are so small that the spindown takes bil- 
lions of years in which case the instability region would 
not really present a no-go area. The second point re- 
quires a thorough understanding of the dynamical evo- 
lution of compact stars and a step towards this goal is 
taken in a companion paper [16 where we show that the 
large amplitude behavior of the bulk viscosity can satu- 
rate r-modes at amplitudes that are large enough for a 
fast spindown. 



^ A potential problem with such an explanation of the data would 
be the observation of a very fast and at the same time very cold 
pulsar, since, after the accretion stops, stars within the stability 
window would re-enter the instability region due to cooling and 
correspondingly would have to spin down to a fraction of their 
Kepler frequency. 
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Figure 7: Instability regions for the different star models considered in this work. Shown are neutron star models with APR 
equation of state (solid), hybrid stars with a small (dot-dashed), a medium (dot-dot-dashed) and a large quark matter core 
(dashed) and strange star models (dotted) with an ideal gas equation of state. The blobs show again the analytic estimates 
for the extrema. Left panel: 1.4 Mq stars. Right panel: 2Mq stars, for which for the considered equations of state stars with 
smaller quark cores could not be found. 
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Appendix A: Analytic form for the R-mode profile 
of constant density stars 



In this appendix we give an analytic solution for the 
density fluctuation profile of the r-mode to leading or- 
der in for a star of constant density. The part of the 
leading order density fluctuation eq. ^ that is generally 
not analytic is given by the change in the gravitational 
potential S^q whose radial part 6^o{r) is deflned by 



TT{m+iy (2m-f-l)! 
It fulflUs the differential equation [5] 
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Here we specialize to the fundamental m — 2 r-mode 
where an analytic solution is possible for the idealized 
case p — const, and dp/dp = const, which is approxi- 
mately realized for strange stars. The solution obtained 
with a computer algebra system reads 



i?3 

— 15L^r) cos( 
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r4((3L2-i?2)sin(f) -3Li?cos(f)) 
in terms of the intrinsic length scale 
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Figure 8: The instability region with direct Urea interactions 
is shown for the 2.21 Mq maximum mass neutron star model 
(solid curve), where direct Urea processes are allowed in the 
inner core but only modified Urea processes in the layer sur- 
rounding it. For comparison the instability region is shown 
when direct Urea reactions are artificially suppressed (dashed 
curve) as well as when they are artificially allowed in the en- 
tire hadronic core (dotted curve). 



Unfortunately the above expression is rather ill behaved 
due to strong cancellations and not suitable for a direct 
evaluation. However, employing the series representa- 
tions of the trigonometric functions it is possible to trans- 
form it into an alternative form 
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Figure 9: Instability regions for the first four multipole r- 
modes (m = 2 to 5) of the different 1.4 Mq star models (top, 
left; neutron star; top, right: hybrid star, small quark core; 
bottom, left: hybrid star, large quark core; bottom right: 
strange star). The minima from eqs. ( |B2[ ) & ( |B3[ ), the second 
minima of the large core hybrid star from eqs. (B8l & (B9| 
and the maxima from eqs. ( B5 1 & ( B6 1 are denoted by the 



dots. For higher multipoles the size of the instability regions 
decreases and they move to higher rotation frequencies, so 
that all modes with m > 7 are entirely stable in the physical 
range of frequencies. In particular in the cases of neutron and 
hybrid stars the right boundaries of the instability regions 
of the higher multipoles are very close to the fundamental 
m = 2 mode, so that several modes can easily be excited if 
the evolution significantly enters the fundamental instability 
region. 
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Appendix B: General semi-analytic expressions for 
the boundary of the instability region 



in terms of hypergeometric functions 
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Correspondingly, \g\ < 0.07 over the entire parameter 
range so that in view of the large uncertainties inherent 
in an r-mode analysis the gravitational potential term 
in eq. ^ can be neglected leaving a simple analytic 
expression for the r-mode profile of strange stars with 
approximately constant density. 



In this appendix we give the general semi-analytic ex- 
pressions for the boundary of the instability region that 
are valid for arbitrary multipoles and can be applied to 
stars involving several shells. The only non-analytic in- 
put in these expressions enters via the radial integral con- 
stants J, S, V and W arising in the gravitational time 
scale, the shear viscosity time scale as well as the low 
and high temperature limit of the bulk viscosity time 
scale. In general these constants, which are given for the 
star models discussed in this work and the normalization 
scales Aqcd — 1 GeV and Aew = 100 GeV in table 



IV have to be computed numerically. However, since 
the dependence of many of the following results on these 
constants is surprisingly weak, these expressions are often 
essentially analytic. 
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Figure 10: Comparison of the instability regions in absolute 
frequencies for the different star models considered in this 
work with the two low mass X-ray binaries Aql X-1 (filled 
square) and SAX J1808. 4-3658 (open circle). The horizontal 
bar gives a partial measure for the error within the model 
computations which should be larger due to uncontrolled as- 
sumptions that are not considered in its size. 



On the low temperature side of the instabihty region 
the shear viscosity dominates and neglecting the bulk 
viscosity yields the analytic result for the boundary in 
this region 
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In the case of a hybrid star with two or more dis- 
tinct shells one of them generally dominates and the 
above expression applies using only the contribution S 
of the dominant shell as well as the corresponding ex- 
ponent (7. The dominant contribution to the inverse 
shear viscosity damping time arises from the hadronic 
shell. Since the shear damping time increases monoton- 
ically with temperature whereas the bulk viscosity de- 
creases in this regime monotonically, the minimum of 
the instability region is taken at the temperature where 
Ts — tb- The minimum frequency is then obtained from 
1/tg + 2/rB = and taking into account that at low 
frequency k sa kq it yields 



Above the minimum the bulk viscosity dominates and 
increases until it reaches a maximum. Within a temper- 
ature range in between these two extrema shear viscosity 
can be neglected and the low temperature approxima- 
tion to the bulk viscosity can be applied. In the general 
case the equation I/tb + I/tg ~ for the boundary can- 
not be solved analytically as a function of T due to the 
non-linear second order terms in 17. However, it is possi- 
ble to solve the equation analytically as a function of 
which provides an equally valid parameterization of the 
segment in between the two extrema 



y^o^ fU„.«o, /2^(2m-f3)(m-f2) 



2m-h2 



(m — 1) 



2™ ^.2 



m(m + l)^™"^((2m + l)!!)^ 



VmA 



9-5 
QCD 



(B4) 



Since there is no analytic expression for the bulk vis- 
cosity around the maximum only for the special case of 
constant density stars a result is possible. In this case 
the maximum frequency is 
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where Ri and Ro are the inner and outer radii of the dom- 
inant shell (0 and R for a homogeneous star) and bars 
denote average quantities over the shell in case the den- 
sity is only approximately constant. The corresponding 
maximum temperature is 
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the dominant shear viscosity at the minimum, in which 
case the above expressions hold, or the bulk viscosities of 
the different layers are identical t^j-* — t'^^ which yields 
analogously 



In general the constant density approximation is only 
valid for quark matter and the corresponding results 
in this case are given in the main text. However, for 
neutron stars the corresponding maximum is in general 
above the Kepler frequency anyway and therefore not 
physically interesting. Far above the temperature where 
the maximum is located the high temperature approxi- 
mation to the bulk viscosity is valid and the boundary 
1/tb + 1/tg = condition yields here 
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For homogeneous stars the present approximations cover 
nearly the entire boundary of the instability region. If 
the star has more than one layer, there may be a second 
minimum. Depending on the size of the layers here it can 
be either the bulk viscosity of the second layer that equals 
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and the corresponding temperature is 
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